Causal perturbation theory in general FRW cosmologies II: 
superhorizon behaviour and initial conditions 
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Abstract 

We describe how causality and energy-momentum conservation constrain 
the superhorizon behaviour of perturbation variables in a general FRW 
■ spacetime. The effect of intrinsic curvature upon the horizon scale is dis- 

cussed and 'white noise' results for the power spectra of generic 'causally- 
seeded' superhorizon perturbations are obtained. We present a detailed 
CO | derivation of the superhorizon energy density power spectrum for curved 

universes, and provide elegant mathematical arguments for a curved uni- 
verse analogue of the Bessel function addition theorem. This yields physical 
insight into both the familiar k 4 result for flat space and the K ^ case: 
. curvature effects are seen to be evident from as early as half the curvature 

^ (— i scale. We consider the implications of these results for the physically and 

Q-f numerically significant issue of setting consistent initial conditions. 
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I. INTRODUCTION 



The standard flat universe inflationary plus dark energy paradigm appears to be 
remarkable accord with recent CMB experiments [p]]. Nevertheless, the confrontation 
with observation remains indecisive, not only because of the significant experimental 
uncertainties, but also because good quantitative accuracy has not yet been achieved 
for competing paradigms such as actively sourced theories, which have so far only been 
extensively studied for cosmic defects in flat FRW (K = 0) backgrounds 0-0. More- 
over, even for flat universes, a subsidiary role for defect networks complementing the 
inflationary power spectrum cannot be excluded ||. The possibilty of searching for dis- 
tinct defect signatures in the CMB, such as nonGaussianity [0] or in the polarisation 
spectra || may go some way toward resolving these remaining uncertainties. However, 
in order to have confidence in cosmological parameter estimation with a purely inflation- 
ary paradigm, it will be necessary to constrain the alternative models, including both 
'primordial' and 'active' perturbation sources, as well as the effects of vector and tensor 
modes and K ^ backgrounds. Here the combination of intrinsic curvature and defect 
sources is particularly interesting, motivating a study of the superhorizon behaviour and 
the choice of consistent initial conditions for such perturbations. Of particular interest is 
the possibility of a testable signature of superhorizon fall-off in 'causal' or 'active' models 
in the CMB polarisation power spectrum. 

The Cauchy problem for cosmological perturbations over a K = FRW background 
has been extensively discussed in the literature. Using energy-momentum conservation, 
and assuming that the early universe was purely homogeneous and isotropic, one can 
show [p|-|l2| that the power spectrum of the stresses should be that of white noise (oc k°), 
while that of the density should be suppressed as k 4 , for wavelengths much larger than 
the horizon; here, k is the comoving wavenumber so this the superhorizon condition is 
k <C aH where H = a/a and a is the scale factor. Assuming that matter undergoes ran- 
dom displacements of length d c at each point in space, while satisfying energy-momentum 
conservation, one can estimate that the density power spectrum turns over to k 4 be- 
haviour at k c ~ 2d~ 1 ||. This result may be used (T^H to simplify the calculations 
to determine consistent initial conditions for numerical analyses. One assumes that all 
perturbations of relevance today come from far outside the initial horizon (kri <C 1) so 
that one may argue that the strong k 4 suppression allows us to set all these variables to 
zero initially. 

For the general K ^ FRW cosmologies, these issues remain unaddressed, although 
there do exist several analyses of the effect of curvature on the CMB assuming primordial 
inflationary perturbations described by power law ansatzes. Some examples include: the 
generalisation of the Sachs- Wolfe effect to include curvature effects [pT| ; refs. and 



T^] treat the open universe case, respectively using a gauge invariant and a synchronous 



gauge formalism — see also ref. |T5[ - and ref. [JT6|] investigated the closed universe 



case, again using a gauge invariant approach. However, the literature does not contain 
(to our knowledge) any direct analysis of the superhorizon effects of energy-momentum 
conservation on active and causal seed sources in intrinsically curved FRW spacetimes. 

In this paper we shall therefore describe the manner in which causality and energy- 
momentum conservation constrain the superhorizon behaviour of the perturbation vari- 
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ables in the general FRW spacetime, in curvilinear coordinates. As discussed in our 
previous paper — ref. [|17]] - we may describe the symmetry properties and stress en- 
ergy conservation laws of these spacetimes in two equivalent languages corresponding to 
integral and differential conservation laws. The language of conserved vector densities 



1^ (c.f. Katz et al |TjJ) uses the conformal geometry of the background spacetime to 
naturally construct integral constraint equations, representing strong conservation laws 
defined over some region of a spacelike hypersurface. These reduce to Noether conserva- 
tion laws when constructed using the Killing geometry of the background. 

Traschen jTj| obtained such integral constraints using as a background manifold a 



maximally symmetric de Sitter spacetime in quasi-Minkowskian coordinates, and these 
were related to the vector densities Ij? by Dereulle et al [pOR . However, in order to make 
use of these integral constraints, we need to know their behaviour on the boundary of 
the integration volume. In the case of bounded (or "local" ) perturbations, the boundary 
terms are assumed to vanish, that is, these previous analyses consider perturbations 
whose stress energy vanishes beyond the horizon. This is a stronger assumption than 
required by the physics, since only the two-point correlators between these quantities 
have to vanish in order to satisfy causality. In other words, a physical assumption has 
been made that is perhaps better suited to the analysis of an explosive event rather than 
more general causal perturbations, such as those seeded by topological defects, which 
may persist beyond the horizon (e.g. cosmic strings). Moreover, the choice of background 
relative to which energy-momentum is formally defined is that of a de Sitter universe, 
rather than the FRW background convenient for most simulations. 

Hence, we turn to the equivalent language of energy-momentum pseudo-tensors |T7| 



outlined in §||. Again, the relative simplicity of the flat space case is complicated by 
the introduction of curvature. In § |111 A| we define the two-point correlator, and discuss 
causality, while the transformed correlator (with respect to the eigenfunctions of the 
Laplacian) is described in §[111 B|. The effect of intrinsic curvature upon the horizon scale 



is dealt with in §1 V A|; and in §|IVB| we obtain results for a general FRW spacetime 



analogous to the white noise power spectrum of two arbitrary functions in flat space. 
Then, in §[V[ we consider in more detail the superhorizon behaviour of the pseudo- 
energy. This not only yields the familiar k A behaviour for the isotropic (OO-)mode on 
superhorizon scales for K = 0, but is also well suited to study the K ^ case. We regain 
the flat space results in the local limit and discuss the effects of intrinsic curvature on 
the isotropic mode in a Helmholtz decomposition of the density perturbation associated 
with causal sources, for both open and closed cosmologies. In we comment on the 
implications of these results for the choice of consistent initial conditions. We work in 
the synchronous gauge because of its ubiquity in numerical simulations and the greater 
physical transparency offered by this gauge choice. 



II. ENERGY-MOMENTUM CONSERVATION 

We shall consider metric perturbations about a general FRW spacetime 

ds 2 = a 2 (7^ + h iiu )dx fl dx v , 
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where the comoving background line element in 'conformal-polar' coordinates (r, \i 0> Q) 
is given by 

juvdafdar = -dr 2 + — - \d X 2 + sin^ x ( d ° 2 + sin2 d( f)\ , ( 2 ) 

\K\ L ^ 

with the function sin^ \ depending on the spatial curvature K as 



sin*- X = I X, K = 0, (3) 




We shall adopt the synchronous gauge defined by the choice h 0fl = 0, so that the trace 
is given by h = h % i (with the convention throughout that Greek indices run from to 3 
and Latin from 1 to 3). 

Here, a = a(r) is the scalefactor, for which we can define the conformal Hubble factor 
7i = a/a, with dots denoting derivatives with respect to conformal time r. The evolution 
of these background quantities is described by the Friedmann equations 

na 2 p = 3(H 2 + K), Ka 2 (p + 3p) = -6H , (4) 

where the background tensor T^ u = (p,p,p,p) includes the dark energy of the universe 



(or cosmological constant). In an earlier paper, |]l7 |, we have shown that one may use 
the Bianchi identities to rewrite the Einstein equations G^ v = R^ v — \g^ v R = kT^ 
(with k = 8tiG), to obtain a description of the perturbed energy-momentum in terms of 
a linear tensor 

KT ° Q = a 2 5G° + Hh-Kh, kt\ = 5G , , 
KT'^S&j-Hltii-hS*.] . (5) 



Note that here we have dropped the deltas in the original definition, understanding that 
the r* 1 ,, are the perturbed parts of the energy-momentum pseudo-tensor only We have 
also separated the perturbed energy-momentum tensor into two parts: the first order 
part ST^v incorporates the stress energy of the radiation fluid, baryonic matter, and cold 
dark matter, and a contribution 0^ representing the stress tensor of an evolving defect 
network or some other causal sources. This is assumed to be small (of order 5T^ U ) and 
'stiff', that is, its energy and momenta are conserved independently of the rest of the 
matter and radiation in the universe and to lowest order its evolution is unaffected by 
the metric perturbations h^ u . 

The components (EH may be simply related to conserved vector densities J M obtained 



by defining energy and momenta with respect to a general FRW background [17|], and 
therefore represent a locally valid definition of these quantities. The conservation equa- 
tions for these quantities are 

Ao + T % = , r% + r J ib . = , (6) 
where the bar denotes covariant differentiation with respect to the spatial three metric 

7ij- 
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The pseudo-energy Tq is a particularly useful variable for tracking the growth of large- 
structure which is why we focus on it in this paper. In ref. Jl7| , we showed that pseudo- 
energy is proportional to the growing mode of the density perturbation on superhorizon 
scales whether the universe is dominated by radiation, matter or curvature (see also 
P^lfl). In addition, we briefly discussed the implications of matching conditions at 
an instantaneous phase transition, with particular emphasis on the choice of consistent 
initial conditions. However, a more thorough treatment requires that we consider not 
the perturbations themselves, but rather their correlators. 



III. TWO-POINT CORRELATORS 

A. Two-point correlators and causality 

Causality is usually defined in terms of two-point correlation functions, which contain 
all the information pertinent to Gaussian perturbations ||22|| . If two random functions b^ v 
and d^ v are assumed to be statistically spatially homogeneous and isotropic, then their 
two point correlation function, C ^ va \{j , t , t') = (6 M!/ (x, r) g?o-a( x + y, r')) is a function 
of time and the 'radial' coordinate |y| = \ y only. The expansion of the universe breaks 
time translation invariance ||, so that we may expect dependence on both r and r', 
rather than just on their difference r — r'. As we are constructing a causal theory of 
cosmological perturbations, we note that causality imposes the constraint 

C^axiy, r > r ') = ( & ^( x > r ) ^a(x + y, r')) = 0, for x y > r + r' , (7) 

so that the two-point correlation function has compact support in y. Here we may assume 
without loss of generality that r > r'. The fact that the correlator has compact support 
in turn means that the Fourier transform (in flat space) of the correlator is analytic, as 
are all the integral transforms with which we shall be concerning ourselves. 

In particular, we shall be concerned with the power spectra of various functions, 
defined as the (Fourier or eigenfunction) coefficient of the equal time two point correlator 
of the function with its complex conjugate: 

C^u{y, r, t) = (6^(x, r)6^(x + y, r)) = — J dV x ft^x, r) fc*„(x + y, r) , (8) 

for an arbitrary function b^ v . Implicit in the final equality of is the "ergodic hypoth- 
esis" which states that one may interchange spatial and ensemble averages. In this case 
the assumption of causality may be expressed as 

<W(y, r, r) = 0, for |y| > 2r = 2 X h ■ (9) 

In the remainder of this paper we shall frequently suppress the (equal) time dependence 
of the correlators. 
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B. The Helmholtz decomposed correlator 



For perturbations over a curved FRW background, it is usual to employ the Helmholtz 
decomposition using the linearly independent eigenfunctions of the Laplacian in polar 
coordinates • We expand all perturbation quantities in terms of the eigenfunctions 



Q( m ), which are the scalar (m = 0), vector (m = ±1) and tensor (m = ±2) solutions to 
the Helmholtz equation 

V 2 Q (m) = 7 ii Q[™ ) = -k 2 Q {m *> , (10) 

where the eigentensor has \m\ suppressed indices (equal to the rank of the perturbation). 
The explicit representations of these eigenfunctions are given in Appendix A. They are 
coveniently labeled in terms of a generalised wavenumber q, and its normalised equivalent 



P = l/ vl-^l) related to k via 



q 2 = k 2 + (\m\ + l)K , <{=► p^ = ~k 2 + (\m\ + l)^- , (11) 

\K 



where k = y\K\k and, in flat space, is the Fourier wavenumber. 

The spectra for flat and open universes (K < 0) are continuous and complete for ft > 0. 
For the open universe, this means that the normalised scalar eigenvalues are greater than 
or equal to unity: k > 1. For the closed (K > 0) case, the spectrum is discrete because 
of the existence of periodic boundary conditons. For scalar perturbations, we then have 
ft — 3, 4, 5, ... since the ft = 1, 2 modes are pure gauge p3|. These spectra yield a set of 



eigenfunctions that is complete, orthonormal, and sufficient to expand all perturbations 
produced in a finite region. 



The literature (see, for instance, refs. p3|j25|j26[1 ) contains several comments about the 
open universe bound k 2 > 1, which could be misinterpreted as implying that perturba- 
tions can only be generated below the curvature scale. However, the particle horizon 
can exceed the curvature scale (see later) which means that, in principle, causal pertur- 
bations also can be generated on these large lengthscales. We shall now briefly place 
these comments in the context of more general eigenfunction decompositions. If one is 
attempting to describe the stochastic properties of inflationary cosmologies it may be de- 
sirable [25,26] to include the eigenfunctions associated with < k 2 < 1 and even k 2 < 0, 



as these may produce effects in the CMB anisotropy patterns — the Grischuk-Zel'dovich 
effect |p7 |. However, this is only necessary if the assumption that the inflaton field is 



in the vacuum breaks down for some correlation length much larger than the curvature 



scale and the effect apparently is not measurable in multi-stage inflation models |28 
Moreover, these eigenfunctions may be obtained only by analytic continuation and are 
not orthonormal, nor are they linearly independent of the k 2 > 1 modes. For square 
integrable functions — such as those describing the correlators of causal perturbations 
— the eigenfunctions labelled by k 2 > 1 are sufficient. 

The k 2 > 1 and k 2 < 1 modes are often referred to as the 'subcurvature' and 'supercur- 
vature modes', because their apparent characteristic wavelength (oc 1/k) is (respectively) 
smaller or larger than the comoving curvature scale. However, this nomenclature is mis- 
leading: the k 2 > 1 modes can contain information about effects above the curvature 
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scale, in which case they are very different to their flat space counterparts. These sub- 
curvature modes are not therefore confined to subcurvature scales. They are complete 
in that they can be used to describe a perturbation on any lengthscale provided it has 
compact support. Now since the horizon can exceed the curvature scale, this complete- 
ness of the k > 1 modes ensures their sufficiency in describing the correlations of causal 
perturbations on all lengthscales and at all times. 

In an FRW universe, the power spectrum (|8|) of any scalar function of y in terms 
of the scalar eigenfunctions = ^(Xx)^m(^xi 0x), where we have decomposed the 
using spherical harmonics (see ( |A11| )). This amounts to a transform on the radial 
coordinate, and a series expansion on the angular part. In particular we wish to expand 
the power spectrum of the contracted second derivative of the pseudo-stresses r given 
in (H) as 

(r*/(x) r ^(x + y)> = yf E/<W) ^{Xy^eyAy)^ ^ r^V(/3), 



Lm 



where 



t 



'/ r P J 9 )UP) = Jl [dV y ^{Xy)YL{e y Ay){T 3] }\-x) r^ix + y)). 



dV x //(x) fdV y r\^^ + y)^{ Xy )Y; m {e y Ay). 

(12) 



Here, in the last expression, we have interchanged integration with respect to x and y. 
In particular, the isotropic 00-mode is given by 



b k (x + y))oo(/3) = ^p/d^ y $?(x y )(r/(x) r^(x + y)) 

= 7=p/ d ^ V"(x) / dV y r\^^ + y)^% Xy 



(13) 

where we have substituted Y 00 = 1/\/Att. 



IV. THE TRANSFORMED HORIZON SCALE AND WHITE NOISE 

A. The transformed horizon scale 

In order to develop intuition as well as to investigate the curved universe superhorizon 
modes in a Helmholtz decomposition, we need to establish the relationship between f3 
and the comoving horizon scale xh — r - This is analogous to the flat space result 
kr fa 3n/2 for a Fourier-Bessel mode with a wavelength of the order of the horizon. We 
therefore consider the expansion, in terms of eigenfunctions of the Helmholtz equation, 
of a 'top hat' window function that vanishes for % x > xh, and is given by A(x) = A 
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for Xx < Xh-, where we choose the arbitrary constant A = l/y/8. The fact that A(x) 
is defined in terms of the radial coordinate only means that we may set I = in the 
Helmholtz decomposition to obtain 



A(x) 



V2 



1 



Now, using dV^ = ^ l f z^ g dxdVL, g = det(^ 



V2^ 
), and appealing to 



J dV A(x)$°( X ) • 
we have 



(14) 



A)o(/?) 



\K\ 3 / 2 y/p 2 + l 



=®Uxh) 



|iY| 3 / 2 /3 2 +l 
1 

|A-|3/2 /3 2_ 1 



sm/SxHCoehxH _ smh cqs ^ 



- cos 



sin cos Pxh 



K < 

= 
K > 



where we have used the explicit expression (|A5|) for $^ given in Appendix A. 

In Fig. |l| we give the transformed 'top hat' functions for two values of the horizon xh, 
normalised to A 0Q — 1 at /3 = 0. When the horizon is less than or equal to the co moving 
curvature scale (xh — 1) ; we see that the open universe transform displays very similar 
behaviour to the flat space Bessel transform. However, on supercurvature scales (xh > 
1), the open and flat universe transforms exhibit different behaviour. In particular, 
we shall choose the first zero of the 'top hat' transform as a means of describing the 
relationship between the horizon scale xh and the corresponding value of the transformed 
variable (3h (the horizon mode wavenumber). From these we will calculate the value of 

" = PhXh , 

in order to discuss how this relationship is affected by curvature. The meaning of this 
relation is that, for a given xh, those modes with f3\H < ol will have characteristic 
wavelengths larger than the horizon. Of course, for flat space and Cartesian coordinates, 
this relation is A = 2n/k < 2ir/kfj = xh- 

The first zeros of A 00 in ( |15D are given by the solutions to 

tanhxij , K < , 

Xh j K = , (16) 
tanxff , K > , 

which are shown in Fig. ^|. Note from ( |I~5| ) that the zeros of ([16]) will be precisely the 
zeros of assuming xh > (for K < 0) and assuming < xh < ^ for the closed 
(K > 0) case. 

The open and flat universes display phenomenologically similar curves. Certainly, as 
can be seen from Fig. |2|, the curves are approximately parallel for xh less than (or of the 
order of) the curvature scale, and that, on these scales, a is therefore the same (~ 37r/2) 
as for flat space: in the local limit, flx = qy — ► kr when x ^ 1 (y ^ \f\K\) 01 when 
\K\ — > 0. However, for comoving distances much larger than xh ~ 4 the effect of the 
curvature is to effectively shrink a as the curvature induced volume growth becomes 
significant. The corresponding limit is 

a = PhXh — ► TT as xh — ► 00 (K < 0) . 

It is clear from this simple discussion that the characteristic lengthscale of a mode cor- 
responding to the eigenvalue k is given by 1//3. 



= tan /3xh — (3 x 



(15) 



8 



FIGURES 




FIG. 1. The transformed 'top hat' function A00: K = (-), K = -1 (- -). 

For an open universe, on very subhorizon or very superhorizon scales, we may simply 
take (3t ^> 1, or (3r <C 1 respectively, which is directly analogous to the flat space 
(kr ^> 1, kr *C 1) limits. However, we note that we must not neglect the precise 
location of the horizon if we are to describe the effect of the curvature scale upon the 
behaviour of the eigenfunctions. We shall return to this issue later (see §[V]), when we 
observe that curvature effects are significant from as early as r = xh — 0.5. 

At this point a brief comment on the difference between the Hubble radius and the 
particle horizon is relevant. In practice, causal perturbations are usually suppressed 
beyond the Hubble radius because this is roughly the maximum distance over which 
correlations can be established in one Hubble time. However, in principle, they can exist 
out to the particle horizon, so we must consider this distinction carefully (e.g. for cosmic 
strings or, in a more extreme case, for exploding neutrino shells). From the Friedmann 
equations (|j) we see that the comoving Hubble radius = (1 — Vt) 1 ^ 2 < 1 for an open 
universe, where the density parameter Vl = pj p cr a £ind p cr ; t = 3H 2 /na 2 . Clearly, this 
implies that the Hubble radius is always less than the curvature scale. In inflationary 
scenarios the Hubble radius is usually the most important quantity and this bound has 
been used to argue that perturbations with characteristic wavelengths smaller than the 
curvature scale are well modelled by a flat universe (see, for example, ref. ||25|| ). However, 
the open universe particle horizon is obtained from the equation 

sinh* (f ) = £ - 1 . (17) 

so that the proper distance du = clXh ~ 2if -1 for O ~ 1, where if -1 is the physical 
Hubble radius. It is therefore possible for the particle horizon to cross the curvature 
scale: Xh > 1 for Q < 0.75. Moreover, in order for xh > 0.5 today (when curvature 
effects become more important — see §[Vp we require only that Q < 0.94 today. This 
is well within current observational constraints, and serves to motivate a more complete 
treatment of the physics around the curvature scale. 
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FIG. 2. Solutions to (pi): K = (-), K = -1 (- -), if = +1 (o) 



For the closed (if > 0) case we see that for the discrete spectrum (3 = 3,4,..., we 
have a ~ 37r/2, which goes over to the (continuous) flat space result in the limit \ — >• 0. 
We note that the notion of supercurvature and superhorizon perturbations with < a, 
Xh ^ 1 does not have meaning in the closed universe. Indeed, from Fig. g, we observe 
that all the (3 have zeros for values of xh < tt/2, at which point the closed universe has 
its maximal circumference and surface area. 



B. The white noise spectrum 

We may use causality and the superhorizon limit to establish that the correlators be- 
tween pseudo-stress must go as a constant on superhorizon scales @Xy ^ a - Substituting 
the small argument series expansion of sin /3x, in the exact expression ( |A4| ) for $° (x), we 
find that the isotropic (/ = 0) mode for the correlator of two arbitrary scalar functions 
may be written as 



2n + l 



00 (_ 1^/3271 f sinh(xy) 

</(x) (x + y)) 00 (/3) = £ 9 ii^i / d ^ < / W »( x + > M X? ^ = 

(2n + l)! 7 I x Wi 



if < 



if > 



On superhorizon scales, we have < a/x y for a given comoving radius Xy ^ Xh 
(where the results of § [IV A| imply that 7r < a < 3n/2). On sufficiently superhorizon 



scales, we may therefore choose (3 <C 1, and the lowest order contribution to ([TJ) will be 
proportional to 0°. Hence, (/(x) g(x + y))oo(/5) oc provided that the integrals are 
finite, which will be the case if we appeal to causality to set < /(x) g(x + y) >= for 
large values of the radial variable Xy This result confirms those by Avelino ||, Pen et 



al || and Veerarghavan and Stebbins [E| for if = 0. 



Moreover, as the vector and tensor eigenfunctions given in Appendix A are propor- 
tional to the scalar radial functions $L we can conclude that correlators of vectors and 
tensors will be premultiplied by a factor behaving as for the scalars. This is in agreement 
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with the results of Durrer et al for K = 0, and corrects the common, but erroneous, 
conclusion that the flat space anisotropic stress power spectra go as k 4 (see, for example, 



refs. POESUSc 



Note that, if we attempt to find the behaviour with respect to k, we run into difficulties 
in describing the higher order corrections. This is a consequence of the fact that there 
will be a constant (curvature) contribution to the k° term from each term (3 2n in the 



series expansion flTgp of $«(x> 



V. SUPERHORIZON BEHAVIOUR OF THE ISOTROPIC MODE 

Of the several arguments in the literature which use energy-momentum conservation to 
yield a k 4 superhorizon power spectrum for flat universes, that provided by Avelino || is 
the most convincing — see also refs. P, pT| , [T2| . Using Cartesian coordinates, and assuming 
that the volume is large enough to justify the ergodic hypothesis, we may equate spatial 
and ensemble averages, and Fourier transform the correlator of the stresses, arriving at 



MM) T ki(k,v')\ = V- 1 I d 3 y e**<7fc(x,i7) 7te(x + y,i/)> . (19) 

J\y\<d c 

Here we have implicitly made use of the convolution theorem, and hence the property 
e jk-y _ e -ik-x e ik-(x+y) Q f exponential kernel used in Fourier theory, while causality 
implies that the correlator has compact support in y. Hence, on superhorizon scales 
(|k|d c <C 1) we may neglect the exponential part of the integrand and find that the 
"power spectrum" of the pseudo-stresses V n . = V\Tij\ 2 / (2n) 3 approachs a constant: we 
have a white noise (k°) spectrum. Noting that partial differentiation with respect to x l 
is equivalent to multiplying the associated Fourier variable by ki, we may use energy- 
momentum conservation (1201) to obtain the result: 



v ™ = V ^r~ = w L dvi L dr]2 L dr]i L dT] > x 

/ d 3 y (^(x,^) 7to(x + y,^)) . (20) 
J\y\<d c 

We conclude that, since (t^Tm) exhibits white noise, then the power spectrum of r 00 
should approach a k A spectrum on superhorizon scales, as has been confirmed in nu- 
merical studies P^| . Similarly, the power spectrum of r° is proportional to k 2 , in this 
limit. 

The difficulty in generalising the above analysis to the general FRW universe lies in the 
lack of a useful equivalent to the exponential partial wave function in the general curved 
(K ^ 0) FRW universe (see ref. f27j), as well as in this argument's dependence upon the 
choice of a Cartesian coordinate system. However, we can recast this argument using the 
eigenmodes of the Laplacian to investigate the constraints imposed by energy-momentum 
conservation and causality in both flat and instrinsically curved FRW universes, regaining 
also the familiar k 4 spectrum in the local (or \K\ — > 0) limit. 
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A. Eigenmode splitting in curved spherical coordinates 



In order to generalise this argument we must be able to split our eigenmodes in 
spherical coordinates under spatial translations to obtain an analogue of the convolution 
theorem. For the flat (K = 0) FRW universe and spherical coordinates we have the 
relationship 

oo 

ja(kr y ) = ^](2n+ l)j n (fcr x )j n (/cr x+y )P n (cosa) , (21) 

ra=0 



where a is the angle between — x and x + y and we have used r_ x = r x [36]. In Appendix 
B we generalise this to the isotropic radial modes $°(x y ) associated with general FRW 
spacetimes: 

oo 

3?(Xy) = E(2^+ l)^(Xx)<^(Xx +y )Pn(cosa) . (22) 

n=0 



The result ( p2j ) is valid for all y. We are, however, principally interested in the 
superhorizon limit, under which this relation simplifies greatly. In the superhorizon 
limit, it is easy to show that $jg(x) °^ (PxY- Also, y superhorizon (fixy ^ 1) implies 
that at least one of — x, x + y is superhorizon also, for any given x. Hence we can make 
use of the first order approximation only: 

K(Xy) « ^(XxK(W , (23) 

for y superhorizon. Thus the lowest radial eigen-mode plays a role analogous to that of 
the exponential in Fourier theory. 

While Appendix B provides a formal proof of the full expansion fl2"2"|), here we just 



obtain (p3|). Consider an eigenfunction expansion of $Jg(x u +v) in terms of $°(x u ), after 
translation to u + v, which is given by 



$°(Xu+v)=E «fm(u)^(Xv)^m^ V) . 
I. rn 

<€» = jrJ dVv $ S(^+v)^(xv)^ m (^,0v) , 

where Ni = V^vr J dVv $l(x v )Yj m (0 v , </> v ) ■ We shall demonstrate that ciq (xu) = 

\^4tv^(Xu), so that (for u = — x and v = x + y) eqn ( ]2"5| ) is confirmed, on superhorizon 
scales where this coefficient carries the leading order behaviour. Here, we give only the 
open universe argument, that for the K > cosmologies being similar. Substituting for 
the volume element, for the I = spherical harmonic, and for each of the two radial 
eigenf unctions using ([Ail), yields 



3 V 71 " f°° f n 

%)(Xu+v) = \ K \3/2Q2 N J d * v sinhx v sin(/?Xv) J d9 sin9 



sm((3x 



u+vj 



\K\y 2 (3 2 N Jo AV AV ^ AV/ y sinhxu+v ' 

Changing variables from 9 to £ = Xu+v using ( |B2|) , we obtain the relation sin9d9 

in(/3Xv)sin(/3Xu 
/3sinhxvsinhxi 



sinh£d£/sinhx v sinhx u , and calculate the inner integral to be 2sm ^ x ^ sm( - /3 ^+ v \ Hence, 
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a o(Xv 



20F sm(p Xu ) roo 
K\^pm Psmhxu Jo Xv PXv 



sin/?Xu 
/5sinhxu 



(24) 



where we have used the normalisation N = {\K\ 3 l 2 l3 2 )~ l J °° dx sin 2 /?x m the last equal- 
ity. Observe that, formally, we should evaluate N , and effect the cancellation in (0), 
in terms of a limit of a volume larger than all length scales of interest. 



B. Integration by parts 



Substituting (|23| ) in ([13]), transforming the inner integration variable to x + y, and 
noting that the integrals are over all space so that d\4 +y = dV y , we see that the I = 
mode is given by the double integral in x and x + y: 



(rVW T^(x + y)> 00 (fc) 



dV x t/(x)$°( X 



V^ 2 

x JdV x+y r 4 /(x + y)$°( Xx+y ). (25) 
In order to evaluate the integrals in ( p5j) we make use of a generalised form of Gauss's 



Law 35 



G 3 



G 2 



W a d^ , 



(26) 



where dV = y~9 {3) dx 1 dx 2 dx 3 and dS a = g (3)ab dS a , dS a = (l/2\)e amim2 d 1 x mi d 2 x m2 . 
Here, we may choose the 3-dimensional vector components diX mi such that dx mi appears 
as the i-th component while the rest vanish (for instance, d2X m ' 2 = (0, dx™ 12 , 0)). As we 
are on a constant time hypersurface, the index a runs over 1,2,3 only, as do the m^. 

Now, if we replace W a with fW a in (|26| ) where / is some scalar quantity (so that 
f-,a = f, a ) and apply the chain rule, we obtain the result 



dvfw a ' a = / dv(fw a y 

G 3 JG 3 



dVf' a W a 
( dS a (fW a ) - / dVf' a W a 

JG 2 JG 3 



(27) 



This is the formal expression of spatial integration by parts on a constant time hypersur- 
face, and a similar result holds with the covariant and contravariant indices interchanged. 
We shall exploit this result twice in each of the two integrals in (f25|). 

We shall proceed in some detail for the integral over x. The x + y integral is similar. 



Since the t*j transform as tensors, setting / 
vector, we find that 



$°(Xx) and Wj = r l so that fWj is a 



\j 



dT4 



dK 
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dK 



- dK 



r 



(2? 



where the third equality is obtained as in (|27|) , and the fourth equality is obtained by 
substituting for Wj, and assuming that the surface integral term vanishes. Explicitly, 
the surface integrals vanish for all but the angular integral over the radial surface at 
infinity, by symmetry, and periodicity. This remaining surface integral may be also be 
set to zero on physical grounds by using the conservation equations (|6]) to obtain 



R 



lim $°JR) / dtt 



T 



J\> 



R 



lim 



R- 



*°JR) 



and by noting that beyond the horizon the functions [$^rjfl (although possibly non- 
zero) should not make a net contribution to the correlator. Physically, we may think of 
this as an expression of the fact that we do not expect there to be any net bulk momentum 
transfers across a radial surface on the very largest scales, since t\ is proportional to 
a well-defined momentum JDJ. This is a much less stringent requirement than that 



previously imposed (i.e. that the functions themselves go to zero beyond the horizon 
P|.|l9|l) which is not expected to be true in most physical situations of interest. Our 
physically motivated and much weaker assumption thus allows for a more complete 
analysis of the superhorizon behaviour. 
Repeating the analysis, we find that 



dK $o (Xx y 



dK 



<(x, 



(29) 



where we have once again ignored the surface terms. Explicitly, the only non-trivial 
surface term is again the radial surface at infinity: 



lim ^jft) / <1Q, 

r—>oc OR 



R 



and, this time, we expect there to be no net anisotropic shears across the radial surface 
on the largest scales, so that this surface term averages to zero and may be ignored. 
Hence, using the connection coefficients, as well as the derivative ( |A7| ) and recursion 
relations for the we see that 



dK <^(Xx)r 
where 



dK (t\ + r 3 3 ) + f dK $2(x) [~~k 2 + 2f K 



T 



1 • 



Sk(P, X) = cot 2 K x - /3 cot K cot(Px) 



(30) 



(31) 



and cot k X — cos kx/ sm xX where cos^X — { cos hx> X-, cos x} f° r {K < 0, K = 0, 
K > 0}. 

Proceeding in a similar fashion for the x+y integral, multiplying to obtain four double 
integrals, and then reversing the steps that led to (EST) — %. e. transforming back to y as 
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the integration variable in the inner integral, and appealing to (|23| 
00-mode as 



we may write the 



/ » b p 



< ? | p 9 )oo(/3) 



( [-P + 2/^/3, Xx)] r\(x) [-P + 2/^/3, Xx +y )] r\(x + y)) 00 (/5) 
+ ( \-k 2 + 2f K (/3, xxjl r\(x) [-f K (p, Xx+y )] (r 2 2 + r 3 3 ) (x + y)) 00 (/3) 



+<-/*(/?, Xx) (r 2 2 + r 3 3 ) (x) [-P + 2f K ((3, x x+y )] r\(x + y)) 00 (/3) 
+ WW, Xx) (r 2 2 + r 3 3 ) (x) [-/*(/?, Xx +y )] (r 2 2 + r 3 3 ) (x + y)) 00 (/3) 



(32) 

Each of the four correlators on the right hand side of (|32|) is equal to a double integral 
similar in form (though, clearly not in integrand) to (P5p. As remarked earlier, for y 
superhorizon, at least one of x, x + y is superhorizon also so that we may take the integral 
over that variable to be finite. Effectively, then, we have obtained the transfromed two- 
point correlators of the contracted derivatives of the pseudo-stresses — and hence the 
pseudo-energy — in terms of functions of k times the transformed pseudo-stresses. 



C. Flat space 



For the flat space (K = 0) case, on sufficiently superhorizon scales (kr 1) we may 
take the limit 



K=Q 



kcot(kr) k 2 

— - w — for kr < 1 



to show that the pre-factors {—k 2 + 2/^ =0 ) 2 , —fK=o(—k 2 + 2fx=o), fx=o are functions 
of k only, in the limit kr — > 0. Hence, we may pull these factors out of the (implicit) 
integrals in fl3"2|) and, taking the limit to infinity, we obtain 



J \j 



<r ...» r p , 19 

\ j\i q\p 



)oo(k) = -[(r\ r\) 00 (k) + (r\ (r 2 2 + r\)) 00 (k) 



r\)oo(k) + ((T 2 2 + T 3 3 ) (r 2 2 + r\)) 00 (k) 



(33) 



where the correlators on the right hand side of ( P3"D are known (see § [i V ri| ) to be propor- 
tional to k° in this limit. Hence, appealing to the conservation equations (|6]) and noting 
that time differentiation does not affect /5-dependence in our eigenfunction expansions, 
we regain the familiar k A power spectrum, as the spectrum for the 00-mode 



<r° (x) T° (x + y))oo(k)<xk 4 . 



(34) 
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D. Open universe 



In this case, there are three regimes of interest: supercurvature, sub-curvature, and 
the case in which x is near the curvature scale. On superhorizon, supercurvature scales, 
X 1, Px ^ 1 which implies that (3 <C 1/x "C 1. Hence, we have that cothv- « 1 and 
~ 1, where these approximations are exact in this limit, and about which we have the 
expansion 

f K <o((3,x) = l-l3cot(/3x) = l-- + ^r--=0(l), 

X 3 

(35) 

Substituting d35|) into (|32|) , we see that the 00-mode tends to a constant on supercur- 
vature, and superhorizon scales. Explicitly, making the k ~ 1 approximation valid on 
these superhorizon and supercurvature scales, we obtain 

(t*/ rV>oo03) = oi<r 1 1 A) 00 (/3) + a 2 (r\ (r 2 2 + r 3 3 )) 00 (/3) 

+«3<(r 2 2 + r 3 3 ) r\)oo(/5) + a 4 ((r 2 2 + r 3 3 ) (r 2 2 + r 3 3 )) 00 (/?) , 

for Xx,y > 1, /?Xx,y < 1 , 

where the are constant factors. Recalling the equations of motion (||) for the energy- 
momentum pseudo-tensor, and that the correlator between two arbitrary functions goes 
as /3°, we see that on very supercurvature and very superhorizon scales, the correlators 
behave as: 

(r° r° ) 00 oc 0° . (36) 

Turning now to the other asymptote, we consider the case in which the horizon is very 
much less than the curvature scale. Then, on superhorizon, and very much subcurvature 
scales, we regain the familiar k 4 power spectrum, up to a constant factor. Using the 
series expansions appropriate for x 1 and fix <S 1, one easily finds that 

f K <o = ^ r = j, / 9 X «l,x«l- (37) 

This has the same form as the superhorizon limit of fx=a- However, one should recall 
that, while (3 (for the K < case) and k (for K = 0) have similar spectra, their meanings 
are different in the two spacetimes — see equation (0). In particular, we have, for the 
open universe with x < 1, two cases: (3 < 1 and (3 > 1. These correspond to mode 
wavelengths approximately equal to and smaller than the curvature horizon, respectively. 
In the former case, fx<o ~ 1/3 and k ~ 1 so that we have the result (EM) once more. 
In the latter, Jk<o ~ /3 2 /3, for large {3, in which case we also have that k 2 w /3 2 . Then, 
recalling the normalisation of fc, we have the local (subcurvature) superhorizon limit: 

fK<o(x,P) » |^/^=o(^,r) . 
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Recalling the conservation equations (|6|), this means that 



<t° (x) r° (x + y)) 00 (/3) oc ~k 4 



\K\ 



P 4 



Xx,y < 1, PXx,y < 1 , 



(38) 



where is the flat space Fourier wavenumber. 

Having dealt with both the extreme cases of local flatness and curvature domination, 
there remains the interesting question of the behaviour near the curvature scale. We 
therefore would like to understand how the k A superhorizon behaviour on sub-curvature 
scales becomes (3° = k° on very super-curvature scales, and the ranges within which 
these approximations are valid. 

Away from the asymptotes dealt with above, fx<o is a function of the radial coordi- 
nate, and therefore is integrated over (twice) in equation (0). This means that we may 
only truly understand the effect of the factors in this equation by explicitly calculating 
the integral, which requires that we specify the pseudo-stresses. As a first approximation, 
however, we may proceed to obtain some understanding of the superhorizon behaviour 
of generic causal perturbations by assuming that the corrections are small enough for us 
to make the approximation 



(x + y) > 00 ~ (PI + P2 + P3 + PA) A 



00 



(39) 



where 



PI = (—k 2 + 2f K<0 ) 2 , 
P3 = -(-& + 2f K<0 )f K<0 



P2 
P4 



f 2 



fi«o(-k 2 



2/ 



K<0) 



for some given x, and where we have assumed that the correlators of the t 1 • are well 



modelled by window functions of the form used in §1 V A| . This assumption amounts to 
ignoring long range correlations, since we assume that Xx ~ Xx+y Hence this approxi- 
mation is valid for Xy small, that is, in the very early universe. 

We provide surface plots for these pre-factors in Figures Q-|5|, from which we can clearly 
see how our approximations are valid only for x < 0.5, and break down as we move into 
the large x an d large f3 regime, doing so much earlier than in the K = scenario. This 
is a geometric effect due to the fact that the majority of the volume in a hyperbolic 
space is located near the radial surface, so that a large 'hyper-sphere' approximates a 
'hyper-shell'. This effect is most noticible on supercurvature scales. However, one can 
calculate that the ratio of hyperbolic volume to surface area rapidly diverges from that 
for a spherical geometry, beginning at around x ~ 0.5. On the other hand, the steep 
slope in the prefactor Pi for large (3 and x < 0.4 is due to the appearance of a constant 
(curvature) term in the (repeated) k factors in this prefactor. However, as x increases, so 
does fx<o, first cancelling, and then overriding this effect. Along the f3 axis in each plot 
(i.e. for x ~ 0), we see that all the pre-factors go as fc 4 /9 for K = and as (f3 2 + l) 2 /9 
for K = —1. Along the x ax is (fl ~ 0), all the pre-factors vanish for K = 0. However, 
for the K = — 1 case, this is not the case, reflecting the small constant (curvature) term. 
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FIG. 6. Log (PI + 2P2 + P3)A00 for K = (-) and K = -1 ( ) 

Computing (r° r° ) using the approximation fl39|) we find the curves plotted in Figure 
^|. These show clearly how the k 4 behaviour on super horizon scales is flattened out, for 
small (3, by the intrinsic curvature. We note that this does not mean that we expect the 
real space correlators to similarly go to a constant on these scales because, of course, 
they vanish by causality. Integrating (|39|) over all /3, with the kernel the constant 
term is suppressed by the measure d/x(/3) = /3 2 dj3. Similarly, calculating the variance 

/*oo rl /*oo 

a 2 = / d/i(/3)(r r° ) 00 (/3)= / dpffliA r° )oo(P) + / cW)(t° r° ) 00 (/3) , 
Jo Jo Ji 

and breaking the integral over all f3 into supercurvature and subcurvature parts, we see 
that the second integral is procedurally the same as for the flat space case, for both 
sub- and superhorizon scales, while the first (containing the constant power spectrum 
constribution) is suppressed by the (3 2 factor in the measure. 

E. Closed universe 

For the closed universe we have the prefactors determined by 

fK>o(P, X) = cot 2 * - P cot x cot(/3x) • (40) 

The difference between the open and closed universe cases is quantified by the differ- 
ences between the behaviour of the function ( pEOP and that of fx<o- To begin with, the 
eigenspectrum for (3 is discrete, and the lowest non-gauge mode is (3 = 3 so that the 
superhorizon limit (3xh 1 can only be attained for xh < 1/3. Of course, for xh 3> 1 
we encounter the periodic properties the closed universe: one cannot have modes that 
are both superhorizon and supercurvature. 

If we take x ^ 1 an d Px ^ 1> then we consider the subcurvature yet superhorizon 
case. Here, we note that cotx w 1/% — %/3, wo we have 

/*>o = ^^ = y, /?X«1, (41) 
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and we regain the flat space results in a similar fashion to the open universe case. Hence, 
an analysis similar to that for the open universe shows that we once again find the local 
limit 

(r (x)r° (x + y))cx^ (42) 

where fc 4 = (ft 2 — l) 2 ~ k 4 /\K\ 2 for large (3, with k being the flat space Fourier wavenum- 
ber. 



VI. DISCUSSION 

It would be of great physical interest if the consequences of this superhorizon be- 
haviour were within the range of observational searches because it would test 'active' 
theories for structure formation, and possibly probe the curvature of the universe. The 
cosmic microwave sky may provide the best prospect, however, for the CMB temperature 
spectrum the superhorizon 'fall-off' expected on the very largest scales will be confused 
by cosmic variance. The most likely signature would be through the CMB polarization 
spectrum because it represents a primordial snapshot of the universe at recombination 
which includes about 10 5 different horizon volumes. The lack of spatial correlations on 
scales corresponding to the horizon size at that time (about 2°), and the fall-off in their 
transformed counterparts, should be observationally testable. Previous analysis of this 
possibility f| , however, has only considered scalar modes, so we comment on the possible 
impact of vector and tensor modes here. We leave discussion of the observational imprint 
of this superhorizon behaviour for a future publication f3(| , concentrating instead on the 



consistent implementation of initial conditions for numerical simulations. 



A. SVT decomposition and superhorizon coherence 

The decomposition of a tensor perturbation into its scalar, vector and tensor parts is 



not without its subtleties [pT2| , since the projection operators yielding the scalar, vector 
and tensor parts are not, in general, analytic. The spatial distribution of the separate ge- 
ometric parts may therefore be significantly different from the physical distribution of the 
entire object, so that the decomposed parts of a causal tensor object can manifest acausal 
features. The transformed variables will then in turn exhibit non-analytic behaviour || , 
complicating the analysis of their superhorizon behaviour. However, causality is merely 
hidden, and, since the arguments presented in this paper for the (scalar) energy density 
perturbation are (trivially) phrased in terms of full tensors, they escape this difficulty 

On the other hand, it is common to employ the assumption of coherence, in which 
case one regains analyticity ||, all the decomposed parts of the perturbation are causal, 
and one may obtain the power spectra by simple power counting. This allows us to write 
down the auto-correlators for both the scalar and the vector parts of the pseudo-tensor 
r^, under the influence of energy- momentum conservation. 

If we decompose the pseudo-tensor t^ v as in ( [All] ), we may write the energy- 
momentum equations (]£]) as 
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Try = -r 



k ' 



T. 



(±1) 
IT 



kf, 



-2k[k 2 -2K]- L t, 



IV 



U(±i) 



2k 2 



(K 



3 Tt 



;(±i) 



V 



(43) 



where we have used the first equation to obtain the final equality in the second. Assuming 
coherence, we find that the scalar (ts) and vector (Try, Ty) degrees of freedom have auto- 
correlators that go as 



(\r S \ 2 )(P) OC k 2 (\T IV \ 2 )(P), 
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r IV \ 2 }((3) oc k i ^\r L \^)- T (7t^m + ^(\Tr\W) 

+ 1QK (<|r L | 2 > - jj(M 2 >(/3)) + lQ^(\r T \ 2 m , 



(l^l 2 )^) oc 



(±1)|2 



>(/*) 



(44) 



on superhorizon scales. Moreover, the correlators for tl, Ty, and Trr have white noise 
superhorizon behaviour (refer to § |IV B and [0]) so we see that in the local limit (x << 1, 
/9 >> 1) we regain the flat space results, with the scalar power spectrum oc fc 4 , while 
those of the vectors is proportional to k 2 . In the supercurvature limit, they all tend to 
a constant. Note that the correlators in ([0]) do not diverge as f3 — > 0, since k tends to 
1 in this limit. 

The results (fPJ) are more general than one might suppose upon first inspection, as an 
arbitrary decoherent source may be written as a 'decoherent sum' of coherent sources, 
each of which will then be amenable to the above arguments Moreover, we believe 
that it is plausible to expect approximate coherence on the very largest scales, since the 
eigenfunctions used to diagonalise the 'decoherent sum' are likely to either all tend to 
white noise on these scales, or to be dominated by just one element in the sum — see also 
ref. [[3(J . This notion has been ventured before as an hypothesis f3lj for scaling sources, 



and is supported by both numerical results |32|], and by the observation that, for causal 
sources, modes with wavelengths much larger than the horizon cannot be significantly 
out of phase with one another. We shall consider this point in more detail elsewhere 



B. Matching conditions and initial conditions 

Finally, we wish to propose a set of initial conditions which would find useful applica- 
tion in numerical simulations of causal perturbation models, such as cosmic defects. We 
shall assume, as is usual, that we are dealing primarily with the evolution of transformed 
quantities. By considered an instantaneous phase transition early in the universe as a 
first approximation to a model for cosmic defects "switching on" , and employing a gauge 
in which constant energy and constant time surfaces coincide, matching conditions can 
be used to show that there exists an entire class of objects which are continuous across 
the transition and are related by gauge transformation to our pseudo-tensor components. 
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In a previous paper (ref. |[T^]), we used the notion of compensation together with a par- 



ticular gauge specification to remove this redundancy such that the t s (pseudo-energy) 
and Ty (divergenceless vector) components of our generalized pseudo-tensor have this 
property: 

M± = o, [t^]± = o. 

While we emphasise that this matching depends on a special gauge choice, it is one which 
should be widely applicable in many physical situations, for example, a defect-forming 
transition . 

For a universe which was unperturbed (and hence homogeneous and isotropic) prior 
to the transition, we may then take t$ = = Ty 11 ^ as natural initial conditions. This 
result is consistent on all scales and is especially appropriate for suppressed superhorizon 
modes. The results of §[V] confirm this conclusion for the density perturbation tracked 
by Ts, and the previous section argues that this conclusion may be extended also to the 
'induced vector' mode t iv and the true vector modes Ty 1 . Hence, if one wishes to set 
initial conditions for perturbations initially outside the horizon, one may consistently 
set both the pseudo-energy r° and the pseudo-momentum r° to zero. If one wishes to 
include long wavelength modes that are affected by curvature, then one can proceed in a 
similar fashion. Although the transformed correlators become constant on these scales, 
at early times they are sufficiently outside the horizon for their amplitude to be strongly 
suppressed. Therefore, the matching conditions we propose as initial conditions will also 
be applicable and self-consistent for a curved universe. 
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APPENDIX A: EIGENFUNCTIONS OF THE LAPLACIAN 

The scalar eigenfunctions are given by solutions to ( |i0|) with \m\ = (c.f., with a 
slightly different normalisation, ref. [p~5j| ) : 



(? lm (x) = \Jl^(y)Yr(0A). (Ai) 

Here the angular dependence is contained in the spherical harmonics which are given 



by Y7™(M) = y ^t 1 fi7S)\ P r( cosd ) eimrp ( where p i m is an associated Legendre function) 
while the radial function $f is 
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2 -L-l 
P ? 



2/3 2 sinhx 

ji(Px) = Jiikr) 
i 

2 



coshx) 



2/3 2 sin X 



-l-l 

p _| +/3 ( cos x) 



if < 
if = 

if => 



(A2) 



with iVg = II^ =0 (/? 2 + n 2 ) and Mg = II^ =0 (/? 2 — n 2 ). By Taylor expanding about fix 
small, we may show that $^ ~ ji(kr) in the limit y — > 0, (3 — > oo. 

The basis {Cpim} is orthonormal as well as complete, and is normalised according to 

(A3) 



^Cp'l'm'Cpim — &(P , P)&l'l&m'm , 

where 5{(3'(3) is the delta function with respect to the measure n{j3): 
J d/i(/3) = / °° p 2 dp , 5(p, F) = js5((3 - p>) 



s(P,P') 



for if < , 
for if > . 



The Z = and / = 1 radial modes are given 
a continuous range of if values): 



by the exact formulae (rewritten for 



sin(/?x) 

P 



sinhx 



if < 
if = 



(A4) 



-±- if > 

sinx 



(P 2 + i)~ 2 [coth x-P cot(px)} if < 



P 



l 

P 

(P 



i 

X 

2 -I 



Pcot{Px) 



if = 



(A5) 



[cot X -/3cot(/3x)] if>0 



The radial function $^ may be recursively calculated using the exact formulae (| 
and the relations: 



(P 2 + i 2 y 



(21 - l)(coth X )$^ - [(3 2 + (I - l) 2 Y 



1-2 



(p 2 ~i 2 y 



(21 - lXcotx)^ 1 - [P 2 -(I- l) 2 ] 2 ^~ 2 



if < 
if = 
if > 



(A6) 



while the derivative with respect to the radial coordinate may be found using 



f Z(cothx)^ - [P 2 + (l + l) 2 ]^ if < 



dx' ! 

Using (|A7|) and 



— =\i&a- P® 



l+l 

P 



if = 



[/(cot^^^-^-^ + l) 2 ]^^ 1 if>0 
for $^ and $^ respectively, one may obtain the relations: 



(A7) 



_d_ 

a v 

_a_ 

d X 
_d_ 



sinh 2 x $l(y) 



sinh 2 X ^(x) K < 



if = . 



(A8) 



Vp^ 



sin 2 X ^(x) if > 
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The vector eigenfunctions satisfy (|T0|) for 
dition D 1 q\ ±1 ^ = 0, where the spectrum is defined by f3 2 
component Qx is given by 



1, subject to the divergenceless con- 
(k 2 + 2K)/\K\. The radial 



QX 



(i) 



i[/(/ + i)] 1/2 $^ m 



-\- K < o 

smn x 

i ' K 

sinx 





fT > 



(A9) 



and similarly for m = — 1, while the remaining (angular) components may be obtained by 
substituting ( |A9| ) into flT0|). The vector eigenfunctions satisfy normalisation conditions 



obtained by replacing (%i> m i(pi 



with r* {±1) r { (±1) in 



with Q; ±1) (/3,/,m) = C^ im . 
2, subject to the transverse-traceless 
(D l Q\j' = = l ) and symmetry (Q^ = Qj' 2 ?) constraints. The spectrum is given 
by j3 2 = (k 2 + 3K)/\K\. The radial-radial component is given by 

i 



.(±i) 



The tensor eigenfunctions satisfy (|T0D for 



m\ 



l(2) 



)(2) i> 



(/ + 2)(z + r 

2/3 2 (/? 2 - 



1/2 



f3 x lm 



sinh 2 y 
J_ 

sin 2 x 



K < 
= 
fT > 



(A10) 



and similarly for m = —2. The remaining components may be obtained by substituting 
( |A10|) in ( |10|) and the constraints. Normalisation has again been achieved by a replace- 
ment of the integrand in (|A|) by Cijpv m > Ca - 
using this decomposition we may decompose arbitrary scalar, vector and tensor per- 
turbations. For example, the pseudo-tensor r^ u defined in (|5]) may be written as: 



T 



x,r 



^■(x,r) 



(All) 



d/i(/3)r 5 g (0) , 

'd/i(/3) [2 (r^-Q^ + W 

+ ^ ) Q (2)i J + rr ) g ( - 2)i J , 

where the coefficients on the right hand side have suppressed dependence on both /3 and 



t. Their explicit form is given in ref. [17]. Here n and represent two 'longitudinal' 
and 'transverse' scalar degrees of freedom, associated with the auxiliary tensors Y jQ^ 

and Q j = k~ 2 Q^^^ + -j 1 -Q^> respectively; tiv is the vector mode 'induced' by 

the auxiliary vector = /•-JvV ( ". 



k- l Q\->; and Tfrp ^ £11*6 the tensor modes 'induced' by 



Q 



(±i) 



-(2k)- 



(±i) 



As well as the transform over the 'radial' coordinate 
/3, there is an implicit sum over indices im which label the spherical harmonics encoding 
the angular dependence. 



APPENDIX B: "ADDITION" RESULTS FOR ¥J X ) 



Here we generalise the result (|2l|) to the K ^ and I 7^ cases, and prove that the 
radial eigenfunctions at u and v are related to those at u + v by 
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nix 



(sinh Xu+v' 



oo 

E 



(sinh Xu) n (sinh x 



(cos cr) 



(Bl) 



where b m is independent of Xu, Xv and a, and may be found by comparing coefficients. 
For u = —x and v = x + y and for the I = case, we may appeal to the common 
normalisation of our radial eigenfunctions (for the K = and K ^ universes) to 
obtain (H). 

Foch |37|] proved this result for the special case I = 0, and for a closed universe only. 
Here we follow the method of Gegenbauer |38| to obtain a more general result for the 
l-th radial eigenmode of an open universe. The proof for the closed universe is similar. 

For curved geometries, the radial coordinates at u and v are related to those at u + v 
via 



COS^Xu+v = COS^XuCOS^Xv 



K 
\K 



■sm^xusm^xvcoscx 



(B2) 



where we use sin^ % an d cos^- % as previously defined, but only for the if / cases. 
Here, a is the angle between u and v (in the hyperplane defined by these vectors), and 
where we may (without loss of generality) choose u pointing towards the north pole such 
that a = 6. For the open (closed) universes, equation ( |B2[ ) clearly reduces to the flat 
space result r^ +v = r 2 + r 2 + 2r u r v cos<r, in the local limit and for \K\ — > 0. 
The n-th radial eigenfunction satisfies the equation 



^ h 2 coth x — 

,2 



dx 2 



dx 



k - 



n(n + 
sinh 2 



nix) = ° • 



(B3) 



From the <K we may construct the quantity 



a, 



%, n(x) 

(sinhx) 1 ' (sinhx)™ 



(B4) 



where v — n+1/2 and defined by ty v = ^/sinh X~nix)i expresses a relationship sim- 
ilar to that between the spherical (half-integer) and ordinary (integer) Bessel functions. 
Then, Q u satisfies the equation 



8 2 Q 8Q 

' - + {2v + 1) coth x-t^- + [k 2 + (u- 1/2) (u + 3/2)]a, 



dx 2 



dx 



. 



(B5) 



Next, we use the addition rule for a hyper-sphere (|B2| ) to write ( p5|) in terms of derivatives 
with respect to Xu and a, assuming without loss of generality that Xv is constant. Hence, 
we obtain the equation 



1 



sinh Xv 



da 2 



d 

2v cot a — 
da 



SK + -T—^ + [2V + 1) COthXuTT^ 
OXu ^Xu 

fc 2 + (z/-l/2)(z/ + 3/2) 



+ 



fl„ = . 



(B6) 



We now assume that VL V can be expanded as Vt v = J2m=n ^mCmi C0S °) where B m is 
independent of a, and the C^(coso") is a (Gegenbauer) polynomial of degree m in cos a. 
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Formally, it may be shown |38|] that + 2z/ cot crj^j C^(coscr) is a constant multiple 
of C^(cosa), so that C^(coscr) may be taken to be the coefficient of a m in the series 
expansion of (1 — 2a cos a + a 2 )~ u . 



From the properties of these Gegenbauer functions C^(coscr), j36fl , we thus arrive at 
the following equation for B m , as a function of Xu : 

d2Rm ' (2u + 1) coth Xu |^ + [k 2 + {y - l/2)(i/ + 3/2)]B m - m{ ? V J m) B m = , 



#Xu ,7 \u 1 ~ "' sinh Xu 

(B7) 

which has solutions B m oc ( -J' nhx ^ , and so, appealing to the symmetry between Xu and 



Xv, we have ([BT]) 
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